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Resumen 


Las crecientes son los gastos máximos anuales de un río y su análisis 
de frecuencias son técnicas estadísticas que permiten estimar valores 
asociados con bajas probabilidades de ser excedidos. Tales 
predicciones permiten el diseño y la revisión hidrológica de las obras 
hidráulicas. La etapa fundamental del análisis de frecuencias consiste 
en seleccionar y ajustar, a los datos o muestra disponible, una función 
de distribución de probabilidades para realizar las predicciones. La 
distribución general de valores extremos (GVE) tiene bases teóricas 
que la convierten en un modelo probabilístico adecuado a crecientes, 
niveles máximos del mar, velocidades extremas de viento, etcétera. 
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Debido a su génesis, la GVE procesa un dato por año, el máximo. Esta 
desventaja fue superada con el método de los r eventos máximos 
anuales. Tal procedimiento se describe y aplica al ajuste de la función 
GVE, a través del método de máxima verosimilitud, utilizando al 
algoritmo Complex de múltiples variables restringidas, para la 
maximización numérica. Se integraron cinco registros de crecientes 
con cinco eventos por año, en estaciones hidrométricas de la Región 
Hidrológica No. 10 (Sinaloa), México. Además, se procesaron dos 
registros de niveles máximos del mar con r = 5, tomados de la 
literatura especializada; uno de ellos es no estacionario al presentar 
tendencia. Se analizaron los resultados y se contrastaron las 
predicciones del ajuste de la GVE con r eventos anuales contra las de 
los métodos clásicos de momentos L y LH. Por último, se formulan las 
conclusiones, las cuales destacan la conveniencia del método descrito 
y señalan las ventajas del uso del algoritmo Complex como técnica 
numérica; debido a lo anterior, se recomienda la aplicación sistemática 


del procedimiento descrito. 


Palabras clave: distribución GVE, método de máxima verosimilitud, 
algoritmo Complex, variables de decisión y dependientes, testde Wald- 
Wolfowitz, error estándar de ajuste, covariables, predicciones. 


Abstract 


Floods are the maximum annual flows of a river and their frequency 
analysis are statistical techniques that allow estimating values 
associated with low probabilities of being exceeded. Such predictions 
allow the hydrological design and review of hydraulic works. The 
fundamental stage of frequency analysis is the selection and fit, to the 


available data or sample, of a probability distribution function to make 
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predictions. The generalized extreme value (GEV) distribution has 
theoretical bases that make ita probabilistic model suitable for floods, 
maximum sea levels, extreme wind speeds, etc. Due to its genesis, the 
GEV processes one data per year, the maximum. This disadvantage 
has been overcome with the method of the r annual maximum events. 
Such procedure is described and applied to the adjustment of the GEV 
function, through the maximum likelihood method, using the Complex 
algorithm of multiple restricted variables, for numerical maximization. 
Five records of floods with five events per year were integrated from 
hydrometric stations of the Hydrological Region No. 10 (Sinaloa), 
Mexico. Also, two records of maximum sea levels with r = 5 were 
processed, taken from specialized literature; one of them is non- 
stationary since it shows a trend. Results were analyzed and the 
predictions of the GEV adjustment with r annual events, were 
compared againstthose ofthe classical methods of moments L and LH. 
Finally, conclusions are formulated, which highlight the convenience of 
the described method and point out the advantages of using the 
Complex algorithm as a numerical technique; due to the above, the 
systematic application of the described procedure is recommended. 


Keywords: GEV distribution, maximum likelihood method, complex 
algorithm, decision and dependent variables, Wald-Wolfowitz test, 
standard error of fit, covariates, predictions. 
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Introducción 


En términos simples, el análisis de frecuencias es una técnica 
estadística que permite ajustar un modelo probabilístico o función de 
distribución de probabilidades (FDP) a una muestra de n datos. Su 
propósito se centra en realizar predicciones, más allá del intervalo de 
los datos empleados para estimar los parámetros de ajuste de la FDP 
seleccionada. Porlo general, las predicciones buscadas son mayores al 
más grande de los eventos observados y, por ello, el análisis de 
frecuencias consiste en el estudio de los eventos pasados con la 
finalidad de establecer probabilidades para ocurrencias futuras. Por 
ejemplo, se estudian los gastos máximos anuales de un río para 
obtener las crecientes de diseño, que son gastos máximos asociados 
con bajas probabilidades de ser excedidos. En resumen, los métodos 
del análisis de frecuencias no predicen el futuro con certeza, pero 
generan modelos probabilísticos que explican y hacen uso eficiente de 
los eventos extremos ocurridos en el pasado (Khaliq, Ouarda, Ondo, 
Gachon, 8: Bobée, 2006). 


El análisis de frecuencias consta de cuatro pasos: (1) verificación 
de la calidad estadística de los datos o muestra disponible. En realidad, 
la validez y exactitud de las predicciones requieren que los datos sean 
independientes y muestren estacionariedad, lo cual implica que hayan 
sido generados por un proceso aleatorio estable en el tiempo; (2) 
selección de una FDP; (3) selección de un método de estimación de los 
parámetros de ajuste de la FDP, y (4) adopción de los resultados o 


predicciones. Comúnmente, se aplica el error estándar de ajuste, que 
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es una medida cuantitativa entre datos observados y estimados con la 
FDP ajustada (Rao 8 Hamed, 2000; Stedinger, 2017). 


En relación con el tópico de la selección de una FDP, la teoría de 
valores extremos es una disciplina estadística que desarrolla modelos 
que describen lo inusual de los datos. Siendo los valores extremos 
escasos por naturaleza, sus predicciones implican fuertes 
extrapolaciones y la teoría de valores extremos proporciona modelos 


que las permiten con base en un fundamento asintótico (Coles, 2001). 


La distribución general de valores extremos (GVE) tiene una base 
teórica que la vuelve imprescindible en el análisis de frecuencias de 
eventos extremos, como crecientes de un río, niveles máximos del 
mar, velocidades extremas de viento, etcétera (Guedes, Soares, 8l 
Scotto, 2004; Khaliq et a/., 2006; An € Pandey, 2007). Sin embargo, 
implícita en su propia formulación de máximo de un lapso o bloque, 
que se expone posteriormente, la distribución GVE emplea un solo 
valor extremo por año y por ello han surgido dos enfoques que buscan 
incorporar más eventos de cada año (Coles, 2001; Khaliq et a/., 2006). 

El primero se denomina de excedencias arriba de un umbral o 
análisis POT (peaks over threshold); utiliza todos los eventos que 
exceden a un valor umbral, por ejemplo, el valor más bajo de los 
máximos anuales observados. Lógicamente, habrá que cumplir con la 
condición de independencia entre los eventos seleccionados. El 
segundo enfoque de uso de más datos se conoce como modelo de r 
estadísticas de orden más grandes (r largest order statistics) y consiste 
en incluir más información de eventos máximos en cada año o bloque 
que cumplan con la condición de independencia entre ellos. Para 
cumplir lo anterior, en la práctica se ha encontrado que r < 5 (Ramesh 
8: Davison, 2002). La ventaja fundamental de ambos enfoques consiste 
en utilizarotros valores máximos de la muestra, en el extremo derecho 
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de la función GVE, para estimar predicciones más precisas o robustas, 
ya que setienen más datos para cuantificarsus parámetros de ajuste 
(Tawn, 1988; Coles, 2001; Khaliq et a/., 2006). 


Los objetivos de este estudio se centran en los cuatro siguientes: 
(1) exponer la teoría en que se basa la aplicación de la distribución 
GVE, con r eventos máximos anuales independientes; (2) describir el 
método de máxima verosimilitud para estimar los tres parámetros de 
ajuste (pu, O, k) de la distribución GVE; (3) integrar cinco registros de 
crecientes con r = 5, dentro de la Región Hidrológica No. 10 (Sinaloa), 
México, para realizar un contraste de predicciones de la distribución 
GVE, entre las obtenidas con el método de r eventos anuales y las de 
los métodos de sextiles, momentos L y momentos LH, estos contrastes 
se realizan con base en el error estándar de ajuste, y (4) procesar con 
el método de los r eventos anuales el registro de niveles máximos del 
mar en Venecia, Italia, que es no estacionario, pues presenta tendencia 
lineal ascendente. 


Teoría operativa 


La distribución GVE 
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En el análisis de frecuencias de datos hidrológicos extremos, como 
crecientes o gastos anuales, niveles del mar, velocidades de viento y 
lluvias máximas, las distribuciones de valores extremos tienen un papel 
preponderante. Su origen se establece a partir de una muestra de 
variables aleatorias (Xi, X2,..., Xn), independientes e idénticamente 
distribuidas (id), que tienen una función de distribución de 
probabilidades F(x). En la práctica, X; son valores del proceso 
observado que se mide cada hora como los niveles del mar o 


diariamente como los gastos de un río. 


El estadístico definido en la Ecuación (1) es el máximo del 
proceso durante n unidades de tiempo de observación (block máxima) 
y en su comportamiento probabilístico se centra el desarrollo de la 


teoría de valores extremos: 


M, = máx[X,, X,,...,X,] (1) 


La distribución de probabilidades acumuladas M, está definida 
por la condición ¡id y es (Smith, 1986; Tawn, 1988; Dupuis, 1997; 
Coles, 2001; An € Pandey, 2007): 


P(M, <x)=P[X,<x,....,X, <x];P(M, < x) =P(X,<x)"""P(X, < 
x); P(M, <x) = [F(0)]" (2) 


Para obtener el comportamiento de la Ecuación (2) conforme n 
tienda a infinito y se eviten soluciones triviales o degenerativas, se 
hace una normalización lineal de la variable M,, con base en las 


constantes an > 0 y bn, la cual es: 
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Me (3) 


Ahora la distribución de los extremos converge a un resultado no 


trivial, definido así: 
P(M; < x) =[F(a, :x+b,]” => G(x) cuando n > co (4) 


La distribución asintótica G(x) debe converger en alguna de las 
tres familias llamadas Gumbel, Fréchet y Weibull. Estas tres 
distribuciones se incluyen en la general de valores extremos (GVE), 


cuya expresión es: 
a -1/x 
0) (- [14 2] ] (5) 


en la cual, y, O y k son los parámetros de ubicación, escala y forma; 
con -o <up<0wo,g>0y-0<k<o, Además, está definida por el 
conjunto (x: [1 + kO<-p4)/0]>0). Si x > 0, la GVE tienetipo Il o Fréchet, 
sin límite superior (uy - 0/k<x< 09). Cuando x < 0, se define la tipo 
III o Weibull, con límite superior (o < x < y - O / k). Finalmente, 
cuando x =0, se llega, en un sentido asintótico, a la distribución tipo I 
o de Gumbel (-o < x < 09), también denominada doble exponencial, 


cuya ecuación es: 


cc = 0-53) 6 
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Seleccionado los máximos (Ecuación (1)) en cada año, las 
predicciones buscadas (X7r) para una probabilidad de excedencia q o 
un periodo de retorno Tren años (Tr = 1/q), se obtienen con base en 
la Ecuación (5), ya que G(x) = 1 - q (Coles, 2001): 


Xp = 4201 [-In(1— q)17") para k + 0 (7) 
Xy, = 4-0" In[-1In(1 — q)] parax=0 (8) 


Conviene aclarar que en el siglo pasado la ecuación de la 
distribución GVE (Smith, 1986; Tawn, 1988; Dupuis, 1997; Hosking 81 
Wallis, 1997; Rao € Hamed, 2000) era: 


G(x) = exp (- [1 —- en (9) 


para la cual, el modelo Fréchet se obtiene con x < 0. Este cambio de x 
por -k en la Ecuación (9) ha sido iniciado, lógicamente por Coles 
(2001), y aceptado por Ramesh y Davison (2002), Katz, Parlange y 
Naveau (2002), y por Khaliq et a/. (2006). Tiene una leve conveniencia 
de signos al manejar la función de verosimilitud y sus restricciones de 
positividad en las ecuaciones (13) a (15) siguientes, y por ello se 
conserva en este estudio. 
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Método de máxima verosimilitud 


De acuerdo con textos consultados (Kite, 1977; Rao € Hamed, 2000; 
Coles, 2001; Kottegoda €: Rosso, 2008; Meylan, Favre, € Musy, 2012) 
el principio de máxima verosimilitud se explica con diferentes niveles 
de detalle. Pero, en general, implica lo siguiente: dada una muestra 
(a, X2,..., Xn) de observaciones ¡¡d, que siguen una distribución Fa con 
parámetros de ajuste 61, 0»,..., Og. Entonces, por definición, la 
probabilidad de obtener un valor X; será: 


P(x;, SX <S x¿+dx¡) = fo(x;): dx; (10) 


siendo fe(xi) la función de densidad de probabilidad. Como los datos 
son jid, la probabilidad de obtener n valores X; será la probabilidad 
conjunta o función de verosimilitud, designada L del inglés /¡kelihood y 


expresada como: 


L(O) = fax)" fal) fon) = M1 fo (xi) (11) 


El método de máxima verosimilitud consiste en encontrar un 
vector Ó de parámetros que hagan máxima a L(60) y por lo tanto a la 
probabilidad de obtener la muestra (X1, X2,..., Xn). Con frecuencia, 
resulta más conveniente tomar logaritmos y trabajar con la función 


logarítmica de verosimilitud (Coles, 2001), es decir: 
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1(6)=1l08L(9) = YE. log f,(x%;) (12) 


lo anterior es aceptable debido a que la función logarítmica es 
monotónica y entonces la función /(0) alcanza su máximo en el mismo 


punto que la función L(0). 


Función /(0) de la distribución GVE 


La función de densidad de probabilidad de la distribución GVE es la 
siguiente (Coles, 2001; Ramesh 8 Davison, 2002): 


f(x; up 0,k) = - [1 +kK ESA - exp (- [1 + K EA (13) 


e + 


donde O > 0, -0 < 4, k< oo y el intervalo de x es tal que [1 + k(x - 
4) / 0] > O, lo cual se indica con el signo + afuera de los paréntesis 
rectangulares. Cualquier combinación de parámetros de ajuste que 
viola la condición anterior de positividad implica que al menos uno de 
los puntos observados (x) está más allá de los puntos finales de la 
distribución, y entonces la función de verosimilitud es cero y la función 
logarítmica de verosimilitud es igual a -oo (Coles, 2001). Si se tiene 
una muestra de máximos anuales Xi, X2,..., Xn, que son 


independientes, la función logarítmica de verosimilitud es: 
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(600 = ola o (14825114 e(208)] 2, [14 (2£0)]% (19) 


Función /(0) de la GVE para r eventos anuales 


Cuando se incluyen otros valores máximos anuales en el ajuste de la 
distribución GVE para buscar una inferencia de predicciones más 
precisa, entonces en cada año se tienen r eventos definidos como: x!? > 
x2>,...,>x”, y lafunción de densidad de probabilidad conjunta permite 
definirla función logarítmica de verosimilitud para n años de r valores 


anuales máximos independientes cada uno (Coles, 2001; Ramesh 81 


Davison, 2002), que es: lu 0) ==n-"-Ino—(1+3) La 2j-aIn|1+ 


a e 


Las restricciones definidas en la Ecuación (13) para los 
parámetros de ajuste se aplican a la expresión anterior y además 
[1 + k(x! — 1)/0] debe ser positivo para cada j = 1,.., r. La Ecuación (14) 
es un caso especial de la Ecuación (15), que se obtiene cuando r = 1. 


Tawn (1988) presenta una ecuación similar a la (15), que considera 
variable cada año el número de eventos r. 
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Función /(0) de la GVE para series no estacionarias 


Los procesos no estacionarios tienen características que cambian 
sistemáticamente a través del tiempo. En los procesos climatológicos 
e hidrológicos, la no estacionariedad se observa por épocas debido a 
los efectos de los patrones climáticos diferentes en varios meses. 
También está presente como tendencia, originada principalmente por 
el cambio climático regional o global (Coles, 2001). 


No existe una teoría general de valores extremos para los 
procesos no estacionarios y, por ello, se sigue un enfoque pragmático 
consistente en utilizar la distribución GVE ajustada por máxima 
verosimilitud, cambiando con el tiempo, o con alguna otra covariable y 
sus parámetros de ajuste (Katz et al., 2002). Para el caso específico 
de una tendencia lineal, se hace variar el parámetro de ubicación con 
el tiempo (Coles, 2001): 


u(t) =PBo+ Bit (16) 


La expresión anterior se sustituye en la Ecuación (15) para 
maximizar tal función y obtener los cuatro parámetros de ajuste 
((Bo PB,,0,k). Otros modelos basados en la GVE para series no 
estacionarias con el enfoque de covariables las describe Coles (2001) 


en su capítulo 6. 
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El algoritmo Complex 


La maximización de la Ecuación (15) para obtener los parámetros 
óptimos de ajuste (y, O, k) de la GVE buscada debe ser abordada de 
manera numérica, y para ello se seleccionó el algoritmo Complex de 
múltiples variables (z) restringidas o acotadas. Su planteamiento 
teórico es el siguiente (Box, 1965): 


Minimizar F(z,,Z,,...,Z,) (17) 


Sujeta a m variables dependientes (y), función de las variables 
de decisión (2): 


Y =D zo): ¿Ya = F(Z Zo 0524) (18) 


Ambas variables tienen límites inferiores y superiores del tipo <, 
es decir, Zing < Z¡ € Zsup Y Yins E Yi E Ysup- El algoritmo Complex es una 
técnica de exploración local, que se guía exclusivamente por lo que 
encuentra a su paso; sus antecedentes, una descripción breve de su 
proceso operativo y su código OPTIM se pueden consultaren Campos- 
Aranda (2003). En Bunday (1985) se tiene otra descripción y código 
de este método de búsqueda. 
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Las designaciones principales en el código OPTIM son NX y NY, 
que definen el número de variables de decisión y dependientes; para 
el caso analizado tres (u, O, k) y n (número de años del registro), 
debido a que las variables dependientes son restricciones de 
positividad, con j = 1,..., r valores por año [1 + k(x!— u)/0]. MI = 500 
es el número máximo de evaluaciones de la función objetivo y NQ = 
25 el número de tales cálculos entre impresión de resultados. Estas 


variables se accesan en la subrutina de lectura de datos. 


Una ventaja importante del código OPTIM radica en permitir un 
fácil acceso de los límites (L = lower, U = upper), nombres y valores 
iniciales de las variables, en la subrutina citada, por medio de las 
designaciones siguientes: XL(D, XU(D, XN$(D, Xx(D, 
YL(J), YU), YNS(0) y YO). Para el caso estudiado, I varía de l a 3, y J 
de 1 a n (número de años del registro). Después, se incluyen los 
criterios de convergencia FA y FR para las deviaciones absoluta y 
relativa de la F. Se utilizaron los valores siguientes: 0.0002 y 0.00001, 


respectivamente. 


La función objetivo se denomina F en el código OPTIM y se accesa 
al final del programa; corresponde lógicamente a la Ecuación (15), con 
nombre FO$="FLMV” de función logarítmica de máxima verosimilitud. 
Se asigna con signo negativo a F debido a que el algoritmo Complex 
minimiza a la función (Ecuación (17)) y se desea maximizara la FLMV. 


Por último, se aclara que por la naturaleza del problema numérico 
planteado, se tiene otra variable, que es NR = 1,..., 5, la cual equivale 
a cada muestra anual de niveles o gastos. La captura de datos se hace 
con dos ciclos anidados: uno para /, variando de 1 a n o número de 
años del registro procesados; y el otro para j, que va de 1 a 5, es decir, 
r de la Ecuación (15). 
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Test de Wald-Wolfowitz 


Esta prueba no paramétrica ha sido utilizada por Bobée y Ashkar 
(1991), Rao y Hamed (2000), y Meylan, Favre y Musy (2012) para 
probar independencia y estacionariedad en registros de gastos 
máximos anuales (xi). Para el caso estudiado, se propuso aplicar la 
prueba al registro de gastos máximos anuales (r = 1), que debe ser 
una muestra de valores aleatorios. Wald y Wolfowitz, basándose en el 
trabajo de Anderson sobre el coeficiente de correlación serial, 


desarrollaron tal prueba, cuyo estadístico es: 


n-1 


R= Y o at 


l= 


(19) 


Cuando el tamaño (n) de la serie o muestra (x;) no es pequeño y 
sus datos son independientes, R procede de una distribución normal 


con media y varianza, dadas por las expresiones siguientes: 


E[R]=R =22% (20) 
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S=SA y SÍ-4SÍ:S¿+4:S1 :S3+5-2S4 52 
Var[R] = =—— == MR (21) 
n-1 (n-1)(n-2) 
en las cuales: 
n 
e y a 
¡=1 
(22) 


Y = rara ei 


El valor de U sigue una distribución normal con media cero y 
varianza unitaria, y se puede usar para probar la independencia de los 
datos de la serie con un nivel de significancia a, por lo común del 5 %. 


En una prueba de dos colas, la variable normal estandarizada es Za/2= 


1.96; entonces, cuando U < 1.96, la serie estará integrada por valores 


independientes (muestra aleatoria). 


Error estándar de ajuste 
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Es el indicador más común (Chai €: Draxler, 2014) para el contraste de 
una distribución de probabilidades a datos reales; se estableció a 
mediados de la década de 1970 (Kite, 1977) y se ha aplicado en México 
haciendo uso de la fórmula empírica de Weibull (Benson, 1962). Ahora 
se aplicará utilizando la fórmula de Cunnane (Ecuación (25)), que de 
acuerdo con Stedinger (2017) conduce a probabilidades de no 
excedencia (p) aproximadamente insesgadas con las distribuciones 
utilizadas en hidrología. La expresión del error estándar de ajuste 
(EEA) es: 


— [E -2p9? ye 
di en 


X¡ son los niveles y gastos máximos anuales (r = 1) ordenados 
de menor a mayor, cuyo número es n y X, los gastos máximos 
estimados con la Ecuación (7) y la probabilidad evaluada con la 
Ecuación (25); np = 3 es el número de parámetros de ajuste de la 
GVE: 


_ 1-0.40 (25) 


— n+0.20 


Búsqueda de la solución óptima y robusta 
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Conforme se avanza de r= 1ar= 5, al aplicar la Ecuación (15) y 
buscar su máximo con el algoritmo Complex, se van procesando cada 
vez más datos, primero n y por último 5n. Esto modifica los valores 
óptimos de los parámetros de ajuste (u, O, k) y también hace cambiar 
al EEA, que se calcula exclusivamente con la serie anual de niveles o 


gastos máximos anuales (r = 1). 


Lo anterior, para poder establecer una comparación objetiva en 
el uso de más información de cada año procesado (solución robusta) y 
buscar una combinación de y, oO y k (solución óptima) que conduzcan 
a un EEA menor o similaral mínimo obtenido con los métodos de ajuste 
por sextiles (Clarke, 1973; Campos-Aranda, 2001); momentos L 
(Hosking 8 Wallis, 1997; Stedinger, 2017), y momentos LH (Wang, 
1997a; Wang, 1997b; Campos-Aranda, 2016). 


Datos procesados 


Registro de niveles en el Mar del Norte 


Guedes Soares y Scotto (2004) presentan una gráfica de cinco niveles 
máximos anuales (metros) en el Mar del Norte, en el periodo 1976- 
1999; por lo cual, n = 24. Sus valores aproximados se muestran en la 
Tabla 1. Estos autores, para asegurarse de utilizar r valores anuales 
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independientes, usan el concepto de duración o longitud estándar de 
tormenta (Tawn, 1988), definida en 480 horas (20 días), al muestrear 


un registro de niveles máximos tomado cada tres horas. 


Tabla 1. Cinco niveles máximos anuales independientes (metros) en 
el Mar del Norte (Guedes Soares 8 Scotto, 2004). 


r valores anuales 


1980 
1981 
1982 
1983 


1984 


1986 
1987 
1988 
1989 


1990 
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1993 12.00 10.98 10.40 


1994 8.46 7.98 7.50 720 
1995 11.48 10.98 10.47 9.98 7.98 


1996 10.97 10.00 8.00 2.00 


cad Mist 


1998 8.52 8.40 8.18 8.04 7.88 


Integración de los registros de crecientes 


Se procesaron cinco registros que fueron integrados con base en la 
información disponible en el sistema BANDAS (IMTA, 2002), en el CD 
1, la cual se denomina “Gastos máximos mensuales”, e incluye día, 
hora y lectura de escala de cada uno de los 12 gastos máximos del año 
analizado. Los cinco registros por procesar de la Región Hidrológica No. 
10 (Sinaloa), México, se exponen de mayor a menor área de cuenca. 
El parteaguas y colector principal de las estaciones hidrométricas: 
Huites, Santa Cruz, Jaina, Guamúchil y El Bledal, se muestran en la 
figura 1 de Campos-Aranda (2014), junto con otras 17 cuencas de tal 
región. 

Para su integración, se aceptó que un lapso de 15 días asegura 
condiciones meteorológicas diferentes en la formación de cada 
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creciente dentro de la Región Hidrológica No. 10 (Sinaloa), México 
(Schulz, 1976; Campos-Aranda, 2000). 


No se realizó una verificación de eventos extremos debido a que 
se acepta que la información hidrométrica que contiene el sistema 
BANDAS (IMTA, 2002) ha sido contrastada y depurada contra los datos 


o mediciones de campo. 


Para buscar que los registros se integren con gastos 
independientes, se siguió el siguiente proceso, en cada año (a) se 
selecciona el gasto máximo anual (x?), y se revisa hacia adelante y 
hacia atrás de tal valor si el gasto mensual adyacente tiene un mínimo 
de 15 días de diferencia con la fecha de este primer máximo; el gasto 
mensual que no cumpla con tal plazo se elimina; después se busca el 
siguiente gasto en magnitud (x?) entre los gastos máximos mensuales 
disponibles(11, 10 o 9 restantes) y se sigue el proceso descrito; así se 
avanza hasta ¡integrar los cinco ¿gastos máximos anuales 


independientes (x). 


Registro de crecientes en la hidrométrica Huites 


La estación de aforos Huites en el río Fuerte de la Región Hidrológica 
No. 10 (Sinaloa), México, con clave 10037 y área de cuenca de 26057 
km2, comenzó a operar en septiembre de 1941 y concluyó en diciembre 
de 1992 (n = 51), cuando inició la construcción de la Presa Luis 
Donaldo Colosio. Con el proceso descrito se integró su registro de cinco 


gastos máximos anuales independientes, que se expone en la Tabla 2. 
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Tabla 2. Cinco gastos máximos anuales independientes (m3/s) en la 
estación hidrométrica Huites, de la Región Hidrológica No. 10 
(Sinaloa), México. 


r valores anuales 


1942 
1943 
1944 
1945 


1946 


2 531.0 
14 376.0 
2 580.0 
1 499.2 


1 164.8 


2 037.6 
3 283.0 
1 262.5 
1 250.0 


445.0 


1 868.8 
1 085.0 
1 024.8 
268.4 
427.5 


718.8 


780.3 
416.4 
768.0 
197.4 
151.0 
634.0 


427.6 
414.2 
474.6 
24,3 
62.2 
205.0 


id IN IN O o 


1948 
1949 
1950 
1951 


1952 


3215.0 
10 000.0 
3 229.3 

677.0 


1 266.0 


799,0 
2.297.5 
1 384.0 

587.5 


895.0 


623.2 
895.3 
961.0 
322.2 
355.8 

68.3 


493.5 
942.4 
439.0 

37.6 
238.3 


a/a 


118.8 
826.4 
438.9 
37.4 
223.5 
39.0 


si Id Ii Ii IE 


1954 
1955 
1956 
1957 


954.8 
4780.3 
695.7 
593.0 


540.1 
1 069.9 
93Lil 


489.0 


481.2 
662.0 
493.2 
380.0 


406.6 
561.8 
278.6 
362.0 


171.4 


35.8 


204.8 
156.2 


1958 3 010.0 1 045.2 894.0 608.5 307.1 
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1959 1 908.0 1 831.0 1 345.5 652.0 544.3 


1960 15 000.0 1 046.0 985.2 721.4 140.5 
1961 1 396.3 905,9 831.6 TIL 2 682.0 
1962 1 620.0 912.0 892.8 501.0 374.5 
1963 2 702.0 1 054.0 980.1 969.2 323.0 
1965 1 944.0 1 787.6 491.2 360.0 147.6 
1966 2 420.0 892.2 738.0 688.9 340.0 
1967 2 505.8 1 310.5 1 192.8 462.5 272.6 
1968 1 534.3 1 118.0 1 019.6 703.8 401.6 


1969 1 508.0 736.0 504.0 417.6 300.0 


1970 1 558.0 1 330.0 970.0 328.7 206.0 


1971 2 200.0 1 176.0 592.0 591.2 576.3 


1972 2 225.0 2 040.0 1142.0 1 109.0 732.1 
1973 7 960.0 2 256.5 1 546.1 1 380.0 372.5 
1974 3 790.0 3 315.0 1 120.0 886.7 587.5 


1975 1 095.0 965.5 532.9 213.2 
1976 2 677.0 1 350.2 1 211.0 Ma al 


1977 1 135.0 622.0 266.9 153.5 120.0 
1978 4 790.0 1 750.0 1119.2 756.0 475.4 
1979 6 860.0 1 001.0 820.0 480.0 302.0 
1980 1 496.0 1 475.0 1 197.8 660.2 452.0 


1981 4 828.1 2 448.0 2 280.0 2 052.0 1 000.3 
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1982 2 450.0 1 571.6 1 085.2 624.0 412.8 


1983 327000 1 439.0 1 400.0 1 006.7 893.0 
1984 5 580.0 1 623.0 1 132.0 624.5 369.2 
1985 3 585.0 1 250.0 1 121.4 925.0 307.6 


1986 1 348.8 1.329,3 436.8 434.0 193.7 


1987 1 429.2 679,2 310.4 


1988 1 866.3 1 494.3 132.6 82.7 47.2 


1989 1 868.5 1 413.9 1 230.3 1 164.7 249.9 
1990 11 558.6 3 544.2 970.1 815.6 376.9 
1991 2001 2 370.0 1721,5 1 517.6 1 266.7 


1992 2.025913 1 564.8 1 348.9 787.3 528.5 


Registro de crecientes en la hidrométrica Santa Cruz 


La estación de aforos Santa Cruz en el río San Lorenzo de la Región 
Hidrológica No. 10 (Sinaloa), México, con clave 10040 y área de cuenca 
de 8919 km?, comenzó a operar en mayo de 1943 y concluyó su 
registro continuo en diciembre de 1980 (n = 38). Con el proceso 
descrito se integró su registro de cinco gastos máximos anuales 


independientes, lo cual se muestra en la Tabla 3 y Figura 1. 
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Tabla 3. Cinco gastos máximos anuales independientes (m3/s) en la 
estación hidrométrica Santa Cruz, de la Región Hidrológica No. 10 
(Sinaloa), México. 


rvalores anuales 


INSI MEU FI 
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Tr = 100 años 


Tr = 50 años 


TF =. 25 años 


Tr= 10 años 


Tr= 5 años 


Gastos máximos anuales (r=1), en m?/s. 


Número de años (í) del registro procesado 


Figura 1. Serie de gastos máximos anuales (r = 1) y rectas de 
predicciones estacionarias en la estación hidrométrica Santa Cruz, de 


la Región Hidrológica No. 10 (Sinaloa), México. 


Registro de crecientes en la hidrométrica Jaina 


La estación de aforos Jaina en el río Sinaloa de la Región Hidrológica 
No. 10 (Sinaloa), México, con clave 10036 y área de cuenca de 8179 
km2, comenzó a operar en enero de 1942 y concluyó su registro 
continuo en diciembre de 1997 (n = 56). Con el proceso descrito se 


integró su registro de cinco gastos máximos anuales independientes, 
lo cual se expone en la Tabla 4. 
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Tabla 4. Cinco gastos máximos anuales independientes (m3/s) en la 
estación hidrométrica Jaina, de la Región Hidrológica No. 10 
(Sinaloa), México. 


PIN A 
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Registro de crecientes en la hidrométrica Guamúchil 


La estación de aforos Guamúchil en el río Mocorito de la Región 


Hidrológica No. 10 (Sinaloa), México, con clave 10031 y área de cuenca 
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de 1645 km2?, comenzó a operar en octubre de 1938 y concluyó en 
diciembre de 1971 (n = 33), cuando inició la construcción de la presa 
Eustaquio Buelna. Con el proceso descrito se integró su registro de 


cinco gastos máximos anuales independientes (Tabla 5). 


Tabla 5. Cinco gastos máximos anuales independientes (m3/s) en la 
estación hidrométrica Guamúchil, de la Región Hidrológica No. 10 


(Sinaloa), México. 


r valores anuales 


1939 299.0 299.0 164.4 [A A 


1940 254.5 37.2 34.1 31.9 
1942 445.0 298.0 161.0 83.0 
1943 1 550.0 1 236.4 298.0 92.3 


1945 916.0 336.0 276.0 228.7 


1947 530.0 133.0 15.5 


1948 648.0 548.0 195.9 40.4 2.8 
1949 375.0 145.4 89.2 72.8 34.0 
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1953 1 173.0 261.3 213.0 MARA 


1954 219.0 115.4 101.2 24.9 
1955 3 507.0 189.0 11.4 2.3 
1956 165.0 148.2 76.4 0.4 
1957 526.0 342.0 57.4 8.2 pen 
Id Ia II Mr Mad 
1959 1610.0 374.0 372.8 22.4 
1960 137.0 130.0 94.5 69.3 
1961 524.5 302.0 134.1 28.5 
1962 985.0 524.0 168.0 112.5 


1963 459.5 311.2 211.5 90.7 


1964 390.0 202.1 128,2 ARA 


1965 449.0 382.5 d.L 1.4 


1966 793.9 687.8 3.4 3.2 
1967 TIOS 329.0 105.7 10.6 
1968 200.0 146.2 132.5 24.8 9.3 


1969 312.0 126.6 84.7 30.3 5.3 


1970 520.0 295.0 256.0 
1971 1 045.0 790.0 175.0 


Registro de crecientes en la hidrométrica El Bledal 
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La estación de aforos El Bledal, en el arroyo del mismo nombre de la 
Región Hidrológica No. 10 (Sinaloa), México, con clave 10027 y área 
de cuenca de 371 km?, comenzó a operar en septiembre de 1937 y 
concluyó en diciembre de 1994 (n = 57). Con el proceso descrito se 
integró su registro de cinco gastos máximos anuales independientes, 


como se expone en la Tabla 6. 


Tabla 6. Cinco gastos máximos anuales independientes (m3/s) en la 
estación hidrométrica El Bledal, de la Región Hidrológica No. 10 


(Sinaloa), México. 


r valores anuales 


1939 eE 204.0 16.5 13.0 
1940 170.0 28.0 19.0 8.4 
1941 118.6 17.1 10.1 4.8 
1942 64.6 33.8 6.4 B.3 


1943 157.2 150.5 32.3 


1944 197.0 46.2 1.0 


1945 414.0 253.2 144.0 [AS 
8.4 0.5 0.5 


1947 174.0 


1948 233.0 141.2 106.0 [ARE 
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1950 44.1 40.4 22.8 

1951 227.4 68.0 29.6 26.7 

1952 169.0 15.4 15.4 [RT 
Roc Mica cil Mal Ml 

1955 152,0 79.8 319 0 

1956 91.4 31.3 48.1 1.0 

1957 283.0 30.5 28.3 0.3 


1959 132.4 29.3 2.4 £.1 0.3 


ps 


1961 122.0 e PURS 3.8 La 0.5 


1962 1 000.0 108.5 56.5 34.7 
1963 335.0 123.0 61.3 0.4 
1964 258.0 189.8 59.0 12.8 
1965 91.5 56.8 ee pies 
— sd 
1967 325.6 115.4 92.7 46.5 
1968 1 576.0 190.4 70.6 6.1 2.0 


1969 228.0 42.0 24.4 3] 0.3 


1971 276.0 260.0 191.0 
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A 
1974 256.1 98.0 64.2 46.9 25.0 
1975 490.0 73.4 30.6 Ll 0.4 
1976 59.6 25.0 25.0 18.4 6.1 
1978 255.0 21.2 5.6 Ol 01 
1979 283.0 32.3 10.8 0.2 0.1 
1980 226.0 ZL1.L 42.6 1.3 0.2 
1981 210.0 118.0 79.8 1.4 0.3 
1982 67.0 62.5 18.0 11,2 4.5 


1984 160.0 153.1 105.2 


1985 278.0 240.0 95.5 86.5 3.2 


1986 T23:0 20.0 15.1 1355 1,9 
1987 107.2 92.4 11.2 2.4 0.3 
1988 234.3 102.5 80.6 2.8 
_ a 
1990 459.3 332.4 86.2 18.1 
19891 73.8 58.2 35.2 10.3 


1992 181.1 161.9 77.6 12.7 


1994 278.2 143.3 101.0 a 
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Registro de niveles con tendencia en Venecia 


En la Tabla 7 y Figura 2 se tienen los primeros cinco niveles máximos 
anuales del mar de los 10 que expuso Smith (1986), para Venecia, 
Italia. 


Tabla 7. Registro de cinco niveles máximos anuales (metros) durante 
51 años en Venecia, Italia (Smith, 1986). 


MecupaMauas 
id ¡Rel pes acia paa 0 100 P05% 


1934 |116| 113 91 |1960|145| 126|123| 116| 114 


1935 |115|107|105|101| 93 | 1961 | 122|108 | 104|100| 100 


il idad Ral hd 


1939 1965 |110|108 | 106|102| 101 


5005 e 
dl a il a ad 
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hs o ad ld Ll il a 
la dd ca a ll 
sal ia a cl ll 
ft a só il il il 


a A 
aid Rd sl ld 


1951 | 151|117|114|109|106| 1977|120| 102| 100 


1952 | 116 91 |1978|132|114 lll 105 


Ed ill hd ad hall he ol 
hd da dad a o a ll 
dd dh al 58M 
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"O. 

220 = 

200 1000 años 

180 00 años Pas 
100 años A 

160 


Niveles máximos anuales, en metros. 


140 
10 años / | ' pS 
120 A e | a de 
| N / N 
MA e Y 
Ñ | VV > 

804 Y decai SE 
Mediana (Tr = 2 años) — 

60 


5 10 15 20 25 30 33 40 45 50 


Número de años (i¡) del registro procesado 


Figura 2. Serie de niveles máximos anuales (r = 1) y rectas de 
predicciones con tendencia lineal y covariable el tiempo (t), en 
Venecia, Italia. 


Resultados 


Prueba de aleatoriedad 


Los resultados del test de Wald-Wolfowitz, definido en las ecuaciones 


(19) a (23), para cada uno de seis registros procesados (Tabla 1, Tabla 
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2, Tabla 3, Tabla 4, Tabla 5 y Tabla 6, columnas 2), son los siguientes: 
Mar del Norte U = -0.518; Huites U = -0.113; Santa Cruz U = -0.827; 
Jaina U = 1.510; Guamúchil U = -1.414; El Bledal U = 0.224, y Venecia 
U = 2.641. Por lo anterior, los primeros seis registros de valores 


máximos anuales (r = 1) son aleatorios y el de Venecia no lo es. 
Ajustes a las series anuales (r = 1) 


En la Tabla 8 se han concentrado los resultados de los ajustes de la 
distribución GVE a cada una de las series anuales de valores máximos 
(r = 1), por medio de los métodos de sextiles, momentos L y momentos 
LH. Tales resultados son los parámetros de ajuste (pu, O, k) y el error 
estándar de ajuste (EEA), cuyo valor mínimo alcanzado en cada serie 


se indica entre paréntesis. 


Tabla 8. Error estándar de ajuste (EEA) y parámetros óptimos de las 


series anuales (r = 1), en los seis registros indicados con los tres 


métodos citados. 


a na 


Mar del Norte 0.189 Ma 048 1.139 -0.3075 


Huites 991.2 Misco 724.345 1 099.256 0.4778 


305.3 814.822 507.620 0.3495 
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Guamúchil 177.2 373.026 260.903 0.3363 


a a ” 


Mar del Norte (0.184) [| 034 1.187 -0.3040 
930.0 Mesas 755.541 1 178.274 0.4459 


Santa Cruz 2002 826.631 402.385 0,4544 
Jaina (190.8) 516.752 381.623 0.4353 


Guamúchil 153.4 362.642 250.526 0.3803 


El Bledal 168.181 123.130 0.2982 
A A 3.636 110.994 16.854 -0.0764 


A > 
Mar del Norte 0.186 10.038 o 106 EM 2391 
(L1) 


Huites (L2) AS 3) 1 Mad 511 al [APR 497 ici 3264 


Santa Cruz Bss 3) Mime 090 Becesad 732 Bis 4044 
(13) 


Jaina (L1) BE 6 [A 368 Man 161 Mad 4296 


Guamúchil (L3) Bull 7) Mica 325 asnall 637 e! 3827 
El Bledal (L4) (32.0) 171.392 116.514 0.3222 
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Venecia (L1) (3.587) 111.690 14.241 0.0405 


Guías para la aplicación del algoritmo Complex 


De inicio se buscó definir límites y valores iniciales de las variables de 
ajuste (1, O, K), y de las dependientes o restricciones de positividad 
que funcionaran para las cinco series por procesar (r = 1,..., 5) de cada 
registro. 


Durante las primeras aplicaciones del código OPTIM surgieron 
varias dificultades, una de ellas estuvo asociada con los límites de las 
variables, que son del tipo <. Entonces, para los límites inferiores no 
se debe poner cero sino 0.10 para los parámetros de ubicación (pu) y 
escala (0), y 0.01 para las variables dependientes. 


Respecto a los límites superiores, se encontró conveniente usar 
como mínimo el doble de los valores calculados para los parámetros yu 
y ad de la GVE de las series de gastos máximos anuales (r = 1) por 
medio del método que condujo al FEA mínimo (Tabla 8). Nunca fue 
necesario cambiar los límites superiores para que el código OPTIM 
funcionara. Para las variables dependientes se adoptó un valor superior 
de 100. 


En relación con los límites del parámetro de forma (k), por lo 
general se usaron como inferior -0.50 y como superior 0.50, o bien -1 


y 1. En pocas ocasiones se tuvo que restringir su límite para obtener 
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resultados en el orden de magnitud del EEA del primer registro 
procesado (r = 1). 

Finalmente, respecto a los valores iniciales de los parámetros de 
ajuste (1, O, k), en general funcionaron bien los valores cercanos al 
mejor ajuste de la GVE (Tabla 8) entre sextiles, momentos L y 
momentos LH con la serie de gastos máximos anuales (r = 1). En la 
Tabla 9 se exponen los límites y valores iniciales utilizados para 


procesar las 30 series de valores máximos anuales. 


Tabla 9. Límites y valores iniciales adoptados para los parámetros de 
ajuste (y, O, k) en la aplicación del algoritmo Complex. 


Registro 


Mar del od 50 
Norte 0.25 


¡E il 10 | 4000 | 1 500 3000 | 1 000 1. il Dec] 


Jaina 0.10 | 1 000 500 0.10 -0.50 | 0.50 | 0.45 


El registro de Venecia se procesó sustituyendo a la Ecuación (16) 


en la Ecuación (15), ahora con cuatro parámetros de ajuste. Los límites 
inferior y superior, y el valor inicial fueron los siguientes: 
fB¿(0.10,200,100), f$,(0.10,2,0.75), 0(0.10,30,15) y k(-0.50,0.50,- 
0.05). 
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Al procesar algunas series de ciertos registros se tuvo que 
cambiar un dato inicial para que el algoritmo Complex llegara a una 
solución cercana a la óptima y no a un mínimo local. Por ejemplo, en 
el registro de Mar del Norte, para la serie de r = 3 se usó -0.05 para el 
k inicial. Para el registro de El Bledal se limitó a 0.40 el valor de 


parámetro de forma k. 
Resultados del algoritmo Complex 


En la Tabla 10 se muestran los principales indicadores de la aplicación 
del algoritmo Complex (aC), como son: la función logarítmica de 
máxima verosimilitud (FLMV), calculada con la Ecuación (15); el error 
estándar de ajuste (E£EA), evaluado con la Ecuación (23), y los valores 
óptimos de los parámetros de ajuste (y, O, Kk). 


Tabla 10. Resultados del algoritmo Complex con cada serie anual de 


r gastos máximos en los registros indicados. 


Registro: == ADN rámetros Sa ¡A 
Inicial final [ERRE 


Mar del AEREA] [2547 9 Dil dl 0.184 E 037 1.124 -0.2712 
Norte 


E 108.8 Nc lla 0.225 10.024 0.995 -0.2001 
EEAmín = 


3| 123.4 154 0.211 10.007 1.029 -0.2164 


0.184 


158.4 0.196 9.993 1.071 -0.2300 
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IO E id Ii o a II It 
_— AAA 
A 

3| 1185.8 | 1182.5 971.5 | 1759.61 | 1193.36 0.3971 


Santa Cruz | 1| 300.9 300.6 111 349.5 856.806 | 457.724 0.3632 
EEAmín = 2| 552.8 552.1 132 267.3 819.749 | 463.835 0.4040 
2198 3| 793.5 787.3 122 408.0 830.896 | 496.459 0.3006 
4| 1033.3 | 1009.9 | 116 364.9 815.842 | 518.103 0.3090 
5| 1295.4 128 470.6 860.041 | 539.006 0.2255 


436.1 435.6 | 127 100.4 505.486 | 372.456 0.4894 


2| 788.2 784.3 230.9 550.876 | 377.100 0.4194 


1106.2 | 1098.4 216.9 550.285 | 376.230 0.4266 


4| 1403.2 | 1394.7 | 115 238.3 547.134 | 380.289 0.4131 


1676.0 | 1663.4 134.9 527.104 | 381.679 0.4593 
Guamúchil 243.6 241.9 153.7 370.294 | 254.203 0.3724 


EEAmín = 2| 448.2 444.9 145 138.0 375.509 | 281.706 0.3518 
iS 3| 630.5 625.2 103 152.6 338.938 | 267.647 0.3607 


4| 802.1 787.2 213 119.3 292.374 | 259.547 0.4991 


o 
684.1 676.5 11 33 144.617 | 116.827 0.3473 
local al hd a st Dt Ea 
ll Bald id cl ad Bal 
leia a a ii Ei 
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También se cita en la columna 1 de la Tabla 10 el EEA mínimo 
obtenido (Tabla 8) con la serie anual de máximos anuales (r = 1), con 
alguno de los tres métodos aplicados: sextiles, momentos L y 
momentos LH. Contra tal valor habrá que comparar los alcanzados por 
el algoritmo Complex, mostrados en la columna 6 (solución robusta). 
Los valores del EEA adoptados se muestran sombreados en la Tabla 


10: son los mínimos alcanzados por el algoritmo Complex. 


Los resultados relativos a las cinco series del registro de Venecia, 
Italia, se han concentrado en la Tabla 11, en la cual se observa que el 
mejor ajuste se logró con r = 2 y un valor de EFA cercano al mínimo 
del método de los momentos LH de 3.587 metros (Tabla 8). 


Tabla 11. Resultados del algoritmo Complex con cada serie anual de 
r niveles máximos en el registro de Venecia, Italia, considerando 


tendencia lineal. 


JN E Parámetros de Wa 


224.6 | 216.1 rr 7.444 Ma 902 Suc 562 | 14.490 | -0.0276 


Discusión 
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Contrastes del EEA 


Con base en los resultados de la Tabla 10 se deduce: (1) para el 
registro de Mar del Norte, que el método de máxima verosimilitud con 
r = 1 iguala al EEA mínimo, alcanzado con el método de momentos L; 
(2) en los registros de las estaciones Huites y Santa Cruz se obtiene la 
mejor opción con r = 5 en la primera y con r = 2 en la segunda; en 
Santa Cruz se obtiene un EEA menor que el mínimo (L3) del método 
de momentos LH; (3) en los registros de las estaciones Jaina, 
Guamúchil y El Bledal se definen como mejores opciones las siguientes: 
r=1,r=4yr = 1, que conducen a un EEA menor que el mínimo 
obtenido con los métodos de momentos L en Jaina y de momentos LH 
en las otras dos estaciones. 


Predicciones del algoritmo Complex 


Seleccionadas, en cada uno de los seis registros procesados, las series 
de r que condujeron a los EEA menores, se contrastan las predicciones 
obtenidas y expuestas en la Tabla 12 contra las del método que 
condujo al menor EEA en la Tabla 8, para establecer las que serán 
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adoptadas, por ser más severas o críticas, las cuales se muestran 


sombreadas. 


Tabla 12. Predicciones para los periodos de retorno indicados, con 
cada serie anual de r gastos máximos independientes en los registros 


indicados. 
Registro Periodos de retorno (Tr) en años 


Mar del 


LALALA ol od 


E 
pe 


Norte 


(metros) 


Huites 
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AAA 
El Bledal 
(m/s) 


67 

81 

32 

72 

12 
211 
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Se aclara que las predicciones mostradas para las cinco series de 
Venecia, Italia, corresponden al final del periodo histórico, es decir, 
utilizando la covariable t con un valor de 51, en la Ecuación (16), al 
aplicarla en la Ecuación (7). 


Asimismo, debido a la amplitud que muestran (33 < n < 57) los 
cinco registros de crecientes procesados de la Región Hidrológica No. 
10 (Sinaloa), México, sus predicciones confiables pueden abarcar hasta 
el periodo de retorno (Tr) de 100 años; sin embargo, se presentan las 
de Tr = 500 y 1 000 años para observar el grado de dispersión que 
tales predicciones robustas muestran en ambos intervalos de 


recurrencia extremos. 


Para el registro de Santa Cruz, en los periodos de retorno (Tr) 
mayores de 50 años se observan en la Tabla 12 predicciones 
ligeramente más grandes con el método de los momentos LH, que con 
las soluciones robustas del método de los r eventos máximos anuales. 
Lo mismo ocurre en el registro de Venecia, pero exclusivamente en los 


periodos de retomo (Tr) de 500 y 1 000 años. 


Por otra parte, para los registros de Huites (r = 5), Santa Cruz 
(r=2) y Guamúchil (r = 4) se obtuvieron predicciones más severas o 
críticas con el ajuste de la GVE y el método de r eventos anuales, que 
las alcanzadas con los métodos clásicos, que utilizan un valor máximo 


anual. Lo anterior es notable en el registro de Guamúchil. 


Finalmente, en los registros restantes (Mar del Norte, Jaina y El 
Bledal) con r = 1, el método de máxima verosimilitud iguala o reduce 
el EEA alcanzado por los métodos clásicos de momentos L y LH. En la 
estación Jaina es factible adoptar los resultados de la serie con r = 5 si 


se desea una solución robusta con un EEA ligeramente mayor. 
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Conclusiones 


El ajuste de la distribución GVE, con el método de máxima verosimilitud 
y utilizando r eventos máximos anuales, es un procedimiento que 
permite utilizar más información de los valores extremos observados 
y, por lo tanto, realiza predicciones más confiables, ya que se están 
estimando sus parámetros de ajuste con más precisión. En resumen, 


conduce a una solución robusta desde un punto de vista estadístico. 


Las siete aplicaciones numéricas descritas mostraron que la 
maximización de la función logarítmica de verosimilitud mediante el 
algoritmo Complex no presentó dificultades computacionales. Además, 
encuentra soluciones óptimas que mejoran el ajuste logrado con los 
métodos clásicos de momentos L y LH, que utilizan únicamente el valor 
máximo anual. Este fue el caso en cuatro de los cinco registros de 


crecientes procesados (Tabla 10). 


Con base en la última aplicación numérica, relativa a los cinco 
niveles máximos del mar en Venecia, Italia, se observa la simplicidad 
para procesar registros no estacionarios que presentan tendencia 
lineal. Lo anterior, debido a la facilidad para incrementar variables de 
decisión, cambiarsus límites y/o valores iniciales, y modificar la función 


objetivo en el algoritmo Complex. 


Con respecto a las predicciones, las cuales variaron de un periodo 
de retomo de 5 a 1 000 años, las soluciones alcanzadas con el 
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algoritmo Complex reportaron en general valores más grandes o 
críticos. Lo anterior fue notable en los registros de Guamúchil, Jaina y 
El Bledal. La excepción fue el registro de la estación Santa Cruz en los 
altos periodos de retorno, donde se obtuvieron predicciones 
ligeramente menores (Tabla 12). 


Tomando en cuenta que no implica gran dificultad integrar 
registros de crecientes independientes con cinco valores máximos por 
año y que la implementación del algoritmo Complex es simple, se 
recomienda aplicar sistemáticamente el ajuste de la distribución GVE 


con r eventos anuales en el análisis de frecuencias de crecientes. 
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